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Recently the PAMELA experiment has released its updated anti-proton flux and anti-proton to 
proton flux ratio data up to energies of ~ 200GeV. With no clear excess of cosmic ray anti-protons 
at high energies, one can extend constraints on the production of anti-protons from dark matter. In 
this letter, we consider both the cases of dark matter annihilating and decaying into standard model 
particles that produce significant numbers of anti-protons. We provide two sets of constraints on 
the annihilation cross-sections/decay lifetimes. In the one set of constraints we ignore any source of 
anti-protons other than dark matter, which give the highest allowed cross-sections/inverse lifetimes. 
In the other set we include also anti-protons produced in collisions of cosmic rays with interstellar 
medium nuclei, getting tighter but more realistic constraints on the annihilation cross-sections/decay 
lifetimes. 



I. INTRODUCTION 

Dark matter (DM) composes approximately 85% of the 
matter density of the universe, but its particle physics 
properties in the WIMP case still remain unknown. Re- 
cently new measurements of cosmic rays|lH7j have gen- 
erated new model building and new constraints on dark 
matter properties. At the same period direct detection 
experiments [8l— Il2l] have provided their own set of con- 
straints on the mass and interaction properties of DM 
particles with nucleons |13l - fl9| |. 

The decay or annihilation products of dark matter gen- 
erally yield equal amounts of matter and antimatter. [70j 
In contrast, antimatter cosmic rays (CR) that are pro- 
duced by Standard model interactions of high energy 
protons and nuclei with the Interstellar Medium (ISM) 
have fluxes that are significantly lower than their mat- 
ter CR antiparticlcs at a broad range of energies. Thus 
antimatter to matter ratios such as the positron fraction 



7(e H 



) or the p/p ratio can provide a very valuable 



probe for identifying DM or constraining its properties. 

[201 ] has provided an update to Q of the anti-proton 
to proton flux ratio measurement by the PAMELA [1| 
satellite, including higher statistics and a new higher en- 
ergy data point extending the measurement of the p/p 
flux ratio up to w 200 GeV. This high ener gy r egion has 
not been probed before with such accuracy, [7l| and thus 
provides us with the opportunity of constraining the pro- 
duction rate of p from DM and to probe new parts of the 
parameter space for the hadronic component of the an- 
nihilation/decay products of DM. In [2l|,[22| the authors 
have provided constraints on annihilating DM via modes 
that produce significant fluxes of ps using the PAMELA 
pup to lOOGeV0. 



In this letter we present the new limits on annihi- 
lating DM cross-sections or decaying DM lifetimes, into 
H^ + V(^ _ , ZZ and bb for a region of DM masses m x . 



II. ASSUMPTIONS 

We consider as reference annihilation mode: \X ~~ ^ 
W + W~ and as reference decay mode: \ — > W + W~ , 
and study the region 100 GeV < m x < 10 TeV for anni- 
hilation and 200 GeV < m x < 20 TeV for decay. We use 
PYTHIA [23J to get the p injection spectra into the ISM 
from the decay and hadronization of the boosted W's. 

High energy electrons and protons loose energy via 
inverse Compton scattering on CMB, IR and starlight 
photons, synchrotron radiation for which the energy loss 
rates ^ scale as 7 2 ,[72| (7 being the Lorentz boost) and 
bremsstrahlung for which ^ ~ 7. Thus high energy p 
and p have much lower energy losses compared to equal 
energy e . On the other hand, high energy protons and 
electrons of energy E 3> m p c 2 have identical rigidities 
and thus diffuse in the ISM in the same way. As a re- 
sult, high energy protons observed locally originate from 
a bigger volume than electrons of the same energy in gen- 
eral. Thus the DM profile distribution is of importance. 
We show the induced constraints from the new p/p flux 
ratio and p flux data for a cuspy profile which we take to 
be Einastof^Jl. 



p(r) = po exp 
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Q 






(1) 



and a cored profile which we take to be a cored isothermal 
profile described by: 



i?| + R 2 C 



(2) 



'Electronic address: ilias.cholis@sissa.it 



The parameters values that we assume are a = 0.17 and 
r_2 = 25 kpc, following Merritt et al. [25| and R c = 



2.8 kpc [26|, |27|. We take the value of the local density 
to be po = 0.4 GeVcm -3 [H] and R e = 8.5 kpc. 

For the propagation of CRs in the Galaxy we use the 
publicly available GALPROP [29| with diffusion within a 
zone of half width of 4 kpc. Diffusion, energy loss, and 
diffusion zone parameters are such that predicted spectra 
of the high energy nuclei fluxes are consistent with local 
CR measurements. 

To determine upper limits on DM annihilation cross- 
sections and decaying DM inverse lifetimes, we first ig- 
nore other astrophysical contributions to the p flux, and 
thus extract upper limits by comparing the PAMELA p 
flux to the p flux from DM. 

Yet, we know anti-protons are produced also by high 
energy protons colliding to ISM protons or nuclei, thus it 
would be an oversimplification to ignore any background. 
Furthermore, anti-protons originating from DM give a 
significantly harder spectrum than that observed at Ep > 
10 GeV (with the exception of m x < 500 GcV(l TcV) for 
annihilation(decay)), failing to fit the observed spectrum 
if they are the dominant component at Ep < 100 GeV. 
Thus we extract a second set of constraints assuming 
the high energy CRs contribution to the p flux. Using 
the same propagation model as the one we assumed for 
the DM originated p fluxes we modified the GALPROP 
assumptions for the pp cross-sections into p [30j, that 
describes the secondary anti-protons component which is 
the dominant component on the flux at high p energies. 
In our assumptions like the authors of [301 ] we fit the ratio 
of p yield from pp reactions to the Simon et al. data [3l( 
in order to have agreement up to the highest T p « 100 



GeV available (where Tp is the p kinetic energy), 
fitted ratio is given by: 



<Tv 
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= 0.12 



GeV 
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GALPROP assumes for the same ratio [3C 



^ = 0A2 (Jjl) 

\GeVj 



-1.67 



1.78. 



(4) 



' pi> 



Our assumptions provide a lower ratio of <Js/a P p than 
[I| at Tp > 30 GeV and slightly higher ratio at 5 - 10 
GeV getting a slightly better agreement to the Simon et 
al. high energy data. 

As lower dark matter masses can contribute in the p 
fluxes down to energies of ~ GeV we allow for the back- 
ground p flux normalization to be decreased up to 10% 
from the best fit value in the absence of a DM originating 
component. 



III. RESULTS 

In figure Q] we present for the case of annihilating 
DM XX -> W+W" the allowed within 95%, 99% and 



99.9% C.L. annihilation cross-section ratios (boost fac- 
tors) (a | v |)/(3.0 x 10- 26 cm 3 s- 1 ). We show the re- 
sults for Einasto profile figure Q}t,b and for cored isother- 
malQJijd. Both limits that come in the absence of a back- 
ground [TJt,c and in the presence of a background [T|3,d 
are shown as described in [TTJ We also give as reference 
the LSP Wino s-wave annihilation cross-section[32| . For 
m x > 800 GeV the constraints are imposed from the 
highest energy data points, while for lower masses the 
data points in the region of 20 — 100 GeV impose the 
constraints on the annihilation cross-section, as the de- 
tails of the p spectrum at those energies start having an 
effect on the fit. 

Since the background that we use fits well the p/p data 
points up to 30 GeV, while gives lower p fluxes at higher 
energies than the PAMELA data suggest, the induced 
constraints on (a \ v |) are more sensitive on background 
assumptions for the lower masses m x . For masses m x > 1 
TeV, the presence of a background results in the allowed 
cross-sections being smaller by a factor of 2, while fol- 
lower masses this factor increases up to ~10. 

The difference in assuming an Einasto versus a cored 
isothermal profile on the allowed cross-sections is a factor 
of less than 2 at m x > 1 TeV, and up to a factor of ^3 at 
the lowest masses, with the Einasto profile setting tighter 
constraints due to its cuspy DM distribution. 

In tableUwe present the allowed within 95% C.L. anni- 
hilation cross-section ratios (a | v |)/(3.0 x 10 -26 cto 3 s -1 ) 
for three different masses, m x = 100 GeV, 1.0 TeV 
and 10 TeV for both Einasto and cored isothermal pro- 
files, and three different annihilation modes: XX ~> ZZ, 
XX ~~ ^ W + W~ , XX ~^ bb. We give the upper limits cal- 
culated both including our background and without it 
(in parentheses). We note also that the limits from the 
new data are -when using only PAMELA data to avoid 
systematic uncertainties between different experiments-, 
a factor of 4 (for XX -► W + W~) to 7 (for XX -> zz ) 
lower for m x = 100 GeV; and a factor of ~ 2 lower for 
the higher masses [74j. 

In figure [2] we present for the case of decaying DM 
X -> W+W~ the allowed within 95%, 99% and 99.9% 
C.L. inverse lifetime ratios r _1 /(5.0 x 10~ 27 s -1 ). As 
with the annihilation scenario we show the results for 
Einasto profile figure [2^i,b and for cored isothermal [2fc,d 
coming from both the absence of a background Gk, c and 
the presence of it [2t>,d- The evolution of the upper limits 
to the mass m x is similar as in the case of annihilation 
(changing m x —> 2m x ). Also as expected the upper lim- 
its in the decay lifetimes -1 are less sensitive to the DM 
halo assumptions than in the annihilation case. 

In table |n] we present the allowed within 95% C.L. 



decay inverse lifetime ratios r /(5.0 x 



for 



m x = 200 GeV, 2.0 TeV and 20TeV. As in table U we 
show results for Einasto and cored isothermal profiles, 
for x -*■ zz -, X ->■ W + W~, x ~> bb, with (without) 
background ps. 

In table IIIII we show the effect of assuming a differ- 
ent diffusion zone thickness (2X) in using GALPROP to 



propagate p in the Galaxy. In each case we used a diffu- 
sion parametrization that gives local CR nuclei, B/C and 
sub Fe/Fe fluxes in good agreement with data (see |33j| ) . 
As wider diffusion zones include propagation of p from 
a bigger volume of the DM halo, they set tighter con- 
straints on annihilation cross-sections and decay inverse 
lifetimes. Other anti-proton backgrounds have also been 
proposed [2ll , |34]-|40J providing upper limits on DM anni- 
hilation/decay rates that can be higher by ~ 2 (see [21|). 
or significantly lower as for the case of j40|, where the 
authors, have studied the significance of the late stages 
of the evolution of supernova remnants in producing p 
(see also [4JJ for the description of their model) . 

In setting the constraints, we ignored DM substruc- 
ture that enhances the p production rate. Including DM 
substructure results in tighter limits on the annihilation 
cross-section or decay lifetime. Since p and p travel from 
greater distances compared to e and the boost from 
substructure in the annihilation case is ~ 1 in the inner 
lOkpc but can grow up to ^10 3 at ^TOOkpc from the 
galactic center [42j| the exact factor by which our results 
can vary by including substructure is a more complex 
question which also depends on the general propagation 
properties assumed at distances ~100kpc from the galac- 
tic center. In the case where the dark disk halo compo- 
nent (seen in recent cosmological simulations including 
baryons [43|, [44J) is significant, the allowed annihilation 
cross-section upper limits can be increased up to a factor 
of 6 0. 

Also even though p = 0.4 GeV cm~ 3 that we have 
used for the local DM density is a value in good agree- 
ment with all estimates, current literature [46l |47] sug- 
gests that it can still be a factor of ~ 2 lower or higher. 
That induces an uncertainty of a factor of ~ 4 (~ 2) to 
the limits that we present for annihilation(decay). 

In figure [3] we show the p/p flux ratio data of [201 fit- 
ted by our background. [75j described in HTl and annihi- 
lating DM, XX -> W + W~ with boost factors BF = 28 
for m x = 1.0 TeV and BF = 2200 for m x = 10 TeV, 
where BF = (a | v |)/(3.0 x 10~ 26 cto 3 s -1 ). We used an 
Einasto profile as given in eq. Q] and our reference propa- 
gation assumptions. We also show for the two DM cases 
presented, the projected AMS-02 [48|, [49| uncertainties 
on the p/p flux ratio after 3 years of data, where the p 
geometrical acceptance is taken to be Ap = 0.25m 2 s r up 
to 300 GeV and dropping to 0.04 at 500 GeV (5^113 ■ 
We note also that systematic uncertainties such as a de- 
teriorating energy resolution with increasing energy and 
at the highest p energies proton spill-over from soft scat- 
terings inside the detector can also increase the p and 
p/p ratio, placing a limit on the highest energies the p 
measurement can be considered accurate. With these un- 
certainties taken properly into account the AMS-02 that 
since May 2011 is on board of the International Space 
Station, with its large acceptance can provide valuable 
information to discriminate among various DM scenarios 
or further constrain them, and together with its measure- 
ments of heavier CR nuclei, decrease the uncertainties in 



Mode 


Profile 


m x = 


100 GeV 


m x = 


= 1.0 TeV 


m x = 10 TeV 


ZZ 

w+w~ 

bb 


Ein. 
Ein. 
Ein. 


1.00 
1.03 
1.4 


(16) 
(15) 
(14) 


90 
68 
81 


(230) 
(160) 

(170) 


5800 (8300) 
2500 (5100) 
3800 (7800) 


ZZ 

w+w~ 

bb 


C.I. 
C.I. 
C.I. 


3.2 
3.6 
3.0 


(23) 
(21) 
(22) 


130 
97 
120 


(340) 
(230) 

(240) 


8500 (12000) 
3600 (7400) 
5600 (11000) 



TABLE I: Allowed upper limits of ratio (a \ v |)/(3.0 x 
10~ 26 cm 3 s~ ) within 95% C.L. for annihilation modes XX ~^ 
ZZ, XX —* W + W~, XX —* bb for Einasto (Ein.) and cored 
isothermal (C.I.) profiles. Values in parenthesis are calculated 
ignoring background p contribution. 



Mode 


Profile 


m x = 


200 GeV 


m x = 2.0 TeV 


m x = 


= 20 TeV 


ZZ 

w+w~ 

bb 


Ein. 
Ein. 
Ein. 


0.11 
0.15 
0.15 


(0.71) 
(0.66) 
(0.67) 


0.43 (1.01) 
0.42 (0.99) 
0.42 (1.04) 


1.6 
1.6 
2.5 


(3.5) 
(3.2) 
(5.0) 


ZZ 

W+W~ 

bb 


C.I. 

C.I. 
C.I. 


0.12 
0.15 
0.15 


(0.74) 
(0.68) 
(0.69) 


0.45 (1.06) 
0.43 (1.01) 
0.43 (1.07) 


1.7 
1.6 
2.5 


(3.6) 
(3.3) 
(5.1) 



TABLE II: Allowed upper limits of ratio r _1 /(5.0x lO" 2 ^ -1 
within 95% C.L. for decay modes x —> ZZ, x — > W + W~ 
X —> bb similar to table [I] 

galactic propagationj50j. 



IV. CONCLUSIONS 

Using the new p PAMELA data [20( we have provided 
a set of constraints on DM annihilation cross-section, de- 
cay lifetimes, for cases where the DM annihilates/decays 
to ZZ, W + W~ , bb producing significant fluxes of high 
energy p. The new constraints are in agreement with 
analysis of the older PAMELA data as in [21| . Uncertain- 
ties in galactic propagation and the shape and local nor- 
malization of the DM halo density allow the constraints 
in the annihilation cross-section/decay lifetimes to vary 
typically by factors of 2-5 but is some cases up to ~ 10. 
For DM models that can explain PAMELA positron frac- 
tion 0, [53[ from either a boosted annihilation rate or an 
appropriate decay rate, the new p data still suggest in 
general DM annihilating/decaying preferentially to lep- 
tons than hadrons 76[ . Based on our most conservative 
upper limits analysis a case of 200 GeV LSP wino as sug- 
gested by [54j is already tightly constrained by the new 
data. Yet wino masses slightly larger as still allowed, thus 



models where DM can annihilate via hadronic modes 55 1 
are still not excluded. Including background antiprotons 
though increases the constrained parameter space to be- 
low 300 GeV, at which point models such as that of [54[ 
can not explain the rise of positron fraction; unless other 
sources such as pulsars are included (for a recent analy- 
sis of constraints on DM annihilating models with a wide 
range of assumed propagation conditions see |56|). Dif- 
fuse gamma-rays from the galactic center and galactic 
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FIG. 1: Upper limits on annihilation cross-section ratios {a \ v |)/(3.0 x 10~ cm s _1 ) vs m x for XX ~ ► W + W~ , for Einasto 
(upper row) and cored isothermal profiles (lower row). Left column: ignoring p background contribution, right column: including 
background ps. Solid lines: 95%, dashed lines: 99%, dotted lines: 99.9% C.L.s. The LSP Wino s-wave annihilation cross-section 
is also given [32| for reference. 



Mode Mass 


L = 2 kpc 


L = 4 kpc 


L = 10 kpc 


W + W~ annih. 0.2 
W+W~ annih. 1.0 
W + W~ annih. 10 


24 (84) 
260 (470) 
7600 (15000) 


6.6 (27) 

68 (160) 

2500 (5100) 


3.1 (11) 

27 (63) 

1000 (2100) 


W + W~ decay 0.4 
W+W~ decay 2.0 
W+W~ decay 20 


0.44 (1.5) 
0.72 (2.9) 

2.8 (7.9) 


0.25 (0.84) 

0.42 (1.03) 

1.6 (4.4) 


0.11 (0.38) 
0.19 (0.56) 
0.69 (1.5) 



TABLE III: Allowed within 95% CL. upper limits of an- 
nihilation to W+W~ cross-section ratio (a | v |)/(3.0 x 
10" 26 cm 3 s _1 ) and decay to W+W~ ratio r _1 /(5.0 x 
10 _27 s _1 ) for varying mass (TeV) and diffusion zone half 
width L. 



ridge j57l |58| and from inner parts of the Galaxy mea- 



sured by F ermi-LAT have been shown to set constraints 
on DM annihilating via hadronic modes (see for instance 
IH-ill). Finally the AMS-02 on board the ISS, due to 
its high acceptance will be able to measure the p/p and 
p spectrum up to at least ~ 300 GeV allowing to further 
constrain or detect DM annihilating/decaying via modes 
that produce significant numbers of hadrons. 
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